Introduction to Statstical Theory of Fluid Turbulence 



Mahendra K. Verm^ 
Department of Physics, IIT Kanpur, Kanpur 208016, India 

This is a brief introduction to the statistical theory of fluid turbulence, with emphasis on field- 
theoretic treatment of renormalized viscosity and energy fluxes. 
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I. INTRODUCTION 

Fluid and plasma flows exhibit complex random behaviour at high Reynolds number; this phenomena is called 
turbulence. On the Earth turbulence is observed in atmosphere, channel and rivers flows etc. In the universe, most 
of the astrophysical systems are turbulent. Some of the examples are solar wind, convective zone in stars, galactic 
plasma, accretion disk etc. 

Reynolds number, defined as UL/v (U is the large-scale velocity, L is the large length scale, and v is the kinematic 
viscosity), has to be large (typically 2000 or more) for turbulence to set in. At large Reynolds number, there are 
many active modes which are nonlinearly coupled. These modes show random behaviour along with rich structures 
and long-range correlations. Presence of large number of modes and long-range correlations makes turbulence a very 
difficult problem that remains largely unsolved for more than hundred years. 

Fortunately, random motion and presence of large number of modes make turbulence amenable to statistical analysis. 
Notice that the energy supplied at large-scales (L) gets dissipated at small scales, say Id- Experiments and numerical 
simulations show that the velocity difference it(x + 1) — u(x) has a universal probability density function (pdf) for 
Id *C I <C L. That is, the pdf is independent of experimental conditions, forcing and dissipative mechanisms etc. 
Because of its universal behaviour, the above quantity has been of major interest among physicists for last sixty years. 
Unfortunately, we do not yet know how to derive the form of this pdf from the first principle, but some of the moments 
have been computed analytically. The range of scales I satisfying Id <C I <C L is called inertial range. 

In 1941 Kolmogorov (19, 2£J,[2l| computed an exact expression for the third moment of velocity difference. He showed 
that under vanishing viscosity, third moment for velocity difference for homogeneous, isotropic, incompressible, and 
steady-state fluid turbulence is 

((un(x + i)- U ||(x)) 3 ) = lm 

where || is the parallel component along 1, (.) stands for ensemble average, and II is the energy cascade rate, which 
is also equal to the energy supply rate at large scale L and dissipation rate at the small scale Id- Assuming fractal 
structure for the velocity field, and II to be constant for all Z, we can show that the energy spectrum E(k) is 



E(k) = K Ko Il 2/3 k- 



■5/3 



where Kko is a universal constant, called Kolmogorov's constant, and L^ 1 <C k <C IJ . Numerical simulations and 
experiments verify the above energy spectrum apart from a small deviation called intermittency correction. 

Availability of powerful computers and sophisticated theoretical tools have helped us understand several aspects of 
fluid turbulence. Some of these theories have been motivated by Kolmogorov's theory for fluid turbulence. Note that 
incompressible turbulence is better understood than compressible turbulence. Therefore, our discussion is primarily 
for incompressible plasma. In this paper we focus on the universal statistical properties of fluid turbulence, which are 
valid in the inertial range. In this paper we will review the statistical properties of following quantities: 



1. Inertial-range energy spectrum for fluid turbulence. 
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2. Various energy fluxes in fluid turbulence. 

3. Energy transfers between various wavenumber shells. 

Many analytic calculations in fluid have been done using field-theoretic techniques. Even though these methods 
are plagued with some inconsistencies, many meaningful results have been obtained using them. Here we will discuss 
items 1-3 in greater detail. 

As mentioned above, pdf of velocity difference in fluid turbulence is still unsolved. We know from experiments 
and simulation that pdf is close to Gaussian for small Su, but is nongaussian for large Su. This phenomena is called 
intermittency. Note that various moments called Structure functions are connected to pdf. It can be shown that 
the structure functions are related to the "local energy cascade rate" H(k). Some phenomenological models, notably 
by She and Leveque 0] based on log-Poisson process, have been developed to compute n(&); these models quite 
successfully capture intermittency in both fluid and MHD turbulence. The predictions of these models are in good 
agreement with numerical results. There are interesting recent advances in intermittency in scalar turbulence [? ]. 
However, here we do not discuss intermittency in great detail. 

Numerical simulations have provided many important data and clues for understanding the dynamics of turbulence. 
They have motivated new models, and have verified/rejected existing models. In that sense, they have become another 
type of experiment, hence commonly termed as numerical experiments. Modern computers have made reasonably 
high resolution simulations possible. The highest resolution simulation in fluid turbulence is on 4096 3 grid (e.g., by 
Gotoh Note that simulations are also used heavily for studying fluid flows around aircrafts and vehicles, in 

atmospheres, engineering devices like turbines etc. 

Fluid turbulence has a larger volume of literature. Here we will list only some of the relevant ones. Leslie |3C| . 
McComb [13, H3, EH , Zhou et al. jfi^l, and Smith and Woodruff Q| have reviewed field-theoretic treatment of fluid 
turbulence. The recent books by Frisch ^3 an d Lesieur j2^| cover recent developments and phenomenological theories. 
Also, the review articles by Orszag [39], Kraichnan and Montgomery j^J, and Sreenivasan are quite exhaustive. 

The outline of the paper is as follows: Section ITfl contains definition of various global and spectral quantities along 
with their governing equations. In Section ITTTI we discuss the formalism of "mode-to-mode" energy transfer rates in 
fluid turbulence. Using this formalism, formulas for energy fluxes and shell-to-shell energy transfer rates have been 
derived. Section HVl contains the existing discussion on Kolmogorov's phenomenology for fluid turbulence. Sections V 
and VI contains very brief introduction on experiments and simulations in fluid turbulence. 

phenomenologies which include Kraichnan's 3/2 mode; Kolmogorov-like models of Goldreich and Sridhar. Absolute 
equilibrium theories and Selective decay are also discussed here. In Section we review the observed energy spectra 
of the solar wind. Section PvTl describes Pseudo-spectral method along with the numerical results on energy spectra, 
fluxes, and shell-to-shell energy transfers. In these discussions we verify which of the turbulence phenomenologies 
are in agreement with the solar wind data and numerical results. In Section IVIII we introduce Renormalization- 
group analysis of turbulence, with an emphasis on McComb's procedure. In Section IVlIII we compute various energy 
fluxes and shell-to-shell energy transfers in fluid turbulence using field-theoretic techniques. Here we also review 
eddy-damped quasi-normal Markovian (EDQNM) calculations of turbulence. 

Section llXl contains a brief discussion on intermittency models of fluid turbulence. The last section (X) contains 
concluding remarks. Appendix El contains the definitions of Fourier series and transforms of fields in homogeneous 
turbulence. Appendix B contains the Feynman diagrams for fluid turbulence; these diagrams are used in the field- 
theoretic calculations. In the last two Appendix (C and D) we briefly mention the mode-to-mode energy transfer 
formalism for scalar and Rayleigh-Benard convection. 
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II. GOVERNING EQUATIONS 



A. Equations for Fluid Dynamics 



Navier-Stokes Equation 



^ + (u- V)u = --Vp th + i/V 2 u+yVV-u, (1) 

where pth is called thermodynamic pressure. The law of mass conservation yields the following equation for density 
field p(x) 



f t + V • {pa) = (2) 
Pressure can be computed from p using the equation of state 

P = I(P) (3) 

This completes the basic equations of fluids. Using these equations we can determine the unknowns (u,p,p). Note 
that the number of equations and unknowns are the same. 

On nondimensionalization of the Navier-Stokes equation, the term Vp becomes (dp/dx 1 ) / p x(C s /U) 2 , where C s 
is the sound speed, U is t he ty pical velocity of the flow, x' is the position coordinate normalized with relative to the 
length scale of the system j4y|. C s — ► oo is the incompressible limit, which is widely studied because water, the most 
commonly found fluid on earth, is almost incompressible (Sp/p <0.01) in most practical situations. The other limit 
C s — > or U C s (supersonic) is the fully compressible limit, and it is described by Burgers equation. The energy 
spectrum for both these extreme limits well known. When U jC s <C 1 but nonzero, then we call the fluid to be nearly 
incompressible; Zank and Matthaeus [sSlfiOl have given theories for this limit. The energy and density spectra are 
not well understood for arbitrary U/C s . 

For most part of this paper, we assume the fluid to be incompressible. In most of the terrestrial experiments, the 
speed of water or air is less than the sound speed. Hence, incompressibility is a good assumption that simplifies the 
calculations significantly. The incompressibility approximation can also be interpreted as the limit when volume of a 
fluid parcel will not change along its path, that is, dp/dt = 0. From the continuity equation (J2J, the incompressibility 
condition reduces to 



V • u = (4) 

This is a constraint on the velocity field u. Note that incompressibility does not imply constant density. However, 
for simplicity we take density to be constant and equal to 1. Under this condition, Eqs. QJGJ reduce to 

— + (u-V)u = -Vp + zA7 2 u, (5) 
Vu = 0. (6) 

When we take divergence of the equation Eq. ©, we obtain Poisson's equation 

-V 2 p = V ■ [(u • V)u] . (7) 

Hence, given u fields at any given time, we can evaluate p. Hence p is a dependent variable in the incompressible 
limit. 

The Navier-Stokes equation is nonlinear, and that is the crux of the problem. The viscous (VV 2 u) term dissipates 
the input energy. The ratio of the nonlinear vs. viscous dissipative term is called Reynolds number Re = UL/u, where 
U is the velocity scale, and L is the length scale. For turbulent flows, Reynolds number should be high, typically 
more than 2000 or so. 
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Table I: Global Quantities in MHD 



Quantity 


Symbol 


Definition 


Conserved in MHD? 


Kinetic Energy 


E 


Jdxu 2 /2 


No 


Kinetic Helicity 


H K 


f dx(u ■ u)/2 


No 


Enstrophy 


n 


J dxuj 2 /2 


No 



B. Energy Equations and Conserved Quantities 



In this subsection we derive energy equations for compressible and incompressible fluids. For compressible fluids 
we can construct equations for energy using Eq. QJ. Following Landau [28] we derive the following energy equation 
for the kinetic energy 



d_ 

dt\2 



-pu + pe 



pu 



V • pu + $ 



(8) 



where e is the internal energy function. The first term in the RHS is the energy flux, and the second term is the 
work done by the pressure, which enhances the energy of the system. The third term,$, a complex function of strain 
tensor, is the energy change due to surface forces. 

In the above equations we apply isoentropic and incompressibility conditions. For the incompressible fluids we can 
choose p = 1. Landau j23| showed that under this condition e is a constant. Hence, for incompressible fluid we treat 
(u 2 )/2 as total energy. For ideal incompressible fluid (v = 0) the energy evolution equation is 



By applying Gauss law we find that 



d_ 

dt 



91 
dt 2 



2 v 



P u 



dx = — 



2 U 



P u 



■dS 



(9) 



(10) 



For the boundary condition u n — or periodic boundary condition, the total energy J 1/2 (u 2 ) is conserved. 

There are some more important quantities in fluid turbulence. They are listed in Table [I] Note that u> is the 
vorticity field. By following the same procedure described above, we can show that in addition to energy, Hk is 
conserved in 3D fluids, while fi is conserved in 2D fluids j^Ell- The conserved quantities play very important role 
in turbulence. 

Turbulent flow contains many interacting "modes", and the solution cannot be written in a simple way. A popular 
approach to analyze the turbulent flows is to use statistical tools. We will describe below the application of statistical 
methods to turbulence. 



C. Necessity for Statistical Theory of Turbulence 

In turbulent fluid the field variables are typically random both in space and time. Hence the exact solutions 
given initial and boundary conditions will not be very useful even when they were available (they are not!). However 
statistical averages and probability distribution functions are reproducible in experiments under steady state, and they 
shed important light on the dynamics of turbulence. For this reason many researchers study turbulence statistically. 
The idea is to use the tools of statistical physics for understanding turbulence. Unfortunately, only systems at 
equilibrium or near equilibrium have been understood reasonably well, and a good understanding of nonequilibrium 
systems (turbulence being one of them) is still lacking. 

The statistical description of turbulent flow starts by dividing the field variables into mean and fluctuating parts. 
Then we compute averages of various functions of fluctuating fields. There are three types are averages: ensemble, 
temporal, and spatial averages. Ensemble averages are computed by considering a large number of identical systems 
and taking averages at corresponding instants over all these systems. Clearly, ensemble averaging demands heavily 
in experiments and numerical simulations. So, we resort to temporal and/or spatial averaging. Temporal averages 
are computed by measuring the quantity of interest at a point over a long period and then averaging. Temporal 
averages make sense for steady flows. Spatial averages are computed by measuring the quantity of interest at various 
spatial points at a given time, and then averaging. Clearly, spatial averages are meaningful for homogeneous systems. 
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Steady-state turbulent systems are generally assumed to be ergodic, for which the temporal average is equal to the 
ensemble average (T^ |. 

Navier-Stokes equation, which is really Newton's equation, is invariant under Galilean transformation 

x = x' + V t' (11) 
t = t' (12) 

where V is the velocity of the primed reference frame with relative to the laboratory frame. Clearly, we can eliminate 
mean velocity of the flow by going to the frame whose velocity is the same as mean velocity of the fluid. Throughout 
this paper we will work in this reference frame. 

As discussed above, certain symmetries like homogeneity help us in statistical description. Formally, homogeneity 
indicates that the average properties do not vary with absolute position in a particular direction, but depends only 
on the separation between points. For example, a homogeneous two-point correlation function is 

( Ui (x,t) Uj (x',t)) = Cij(x- x',t) = Cij(r,t). (13) 

Similarly, stationarity or steady-state implies that average properties depend on time difference, not on the absolute 
time. That is, 

(« i (x,t)« i (x,*')) =Gy(x,t-0- (14) 

Another important symmetry is isotropy. A system is said to be isotropic if its average properties are invariant under 
rotation. For isotropic systems 

(« i (x,t)« i (x / ,t))=C ii (|x-x / |,t) = a j -(|r|,t). (15) 

Isotropy reduces the number of independent correlation functions. Batchelor 1] showed that the isotropic two-point 
correlation could be written as 

dj (r) = C7« {r) nr] + {r)5 l3 (16) 

where C« and are even functions of r = |r|. Hence we have reduced the independent correlation functions to 
two. For incompressible flows, there is only one independent correlation function 

In turbulent fluid, fluctuations of all scales exist. Therefore, it is quite convenient to use Fourier basis for the 
representation of turbulent fluid velocity and magnetic field. Note that in recent times another promising basis called 
wavelet is becoming popular. In this paper we focus our attention on Fourier expansion, which is the topic of the 
next subsection. 



D. Turbulence Equations in Spectral Space 

Turbulent fluid velocity u(x, t) is represented in Fourier space as 

/dk 
^-pu(k,t)exp(ik-x) (17) 

u(k,i) = / dxu (x, t) exp (— ik ■ x) (18) 



where d is the space dimensionality. 

In Fourier space, the equation for incompressible fluid is 0] 



d_ 

01 



(19) 



+ vk 2 ^j ui(k, t) = -ikip(k,t) - ikj J j0^itj(k- p,t)iii(p,t) 
with the following constrains 

k-u(k) = 0, (20) 

The substitution of the incompressibility condition k • u (k) =0 into Eq. l(T9jl yields the following expression for the 
pressure field 



/i \ fcj^kj F dip r /i \ / m 

p( k > = — -gr J 7^w[^(k-p,%(p,«)] 



(21) 
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Table II: Various Spectra of Fluid Turbulence 



Quantity 


Symbol 


Derived from 


Symbol for ID 


Kinetic energy spectrum 


C(k) 


MkK(k')) 


E(k) 


Enstrophy spectrum 


fi(k) 


Mk) Wj (k')} 


fi(Jfe) 



Note that the density field has been taken to be a constant, and has been set equal to 1. 

It is also customary to write the evolution equations symmetrically in terms of p and k — p variables. The sym- 
metrized equations are 



d_ 

dt 



+ vk 1 Ui(k,t) 



^|m( k ) J ^[%(p,tK(k^P 



,*)] 



(22) 



where 



^n(k) = %P im (k)+fc m P«(k); 

Pirn (k) = 5fc 



hi km 



k 2 ' 

Energy and other second-order quantities play important roles in fluid turbulence. For a homogeneous system these 
quantities are defined as 

( Ui (k,t) Uj (Vi',t)) = C7 l ,(k,t)(2^ !i <5(k + k'). 
The spectrum is also related to correlation function in real space 

/dk 
^-pCy (k)exp (ik-r). 

For isotropic situations we can take CV,(k) to be an isotropic tensor, and it can be written as [1] 

CM = PMCik). (23) 

When turbulence is isotropic, then a quantity called ID spectrum or reduced spectrum E(k) defined below is very 
useful. 

H<" 2 > -i /<^ c «< k > 

E{k)dk 



1 



where Sa 



, dk f s . P»(k)C(k) 
2 / (27r) d K ' y ' 

2TT d / 2 /T(d/2) is the area of d— dimensional unit sphere. Therefore, 

-d-i Sd(d - 1) 



E(k) = C(k)k° 



2(27r) rf 



(24) 



Note that the above formula is valid only for isotropic turbulence. We have tabulated all the important spectra of 
fluid turbulence in Table ITfl 

The global quantities defined in Table U are related to the spectra defined in Table ITT1 bv Perceval's theorem 0]. 
Since the fields are homogeneous, Fourier integrals are not well defined. In AppendixElwe show that energy spectra 
defined using correlation functions are still meaningful because correlation functions vanish at large distances. We 
consider energy per unit volume, which are finite for homogeneous turbulence. As an example, the kinetic energy per 
unit volume is related to energy spectrum in the following manner: 



i / d*l (u 2 ) = I 
L d 2 x 1 2 



dk 



(27T) 



^(k) = / E(k)dk 
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Similar identities can be written for other fields. 

In three dimensions we have another important quantities called kinetic helicities. In Fourier space kinetic helicity 
iSRr(k) is defined using 

( Ui (k, t) il, (k', t)) = P it (k) i^(k)(2^) d <5(k + k') 
The total kinetic helicity Hk can be written in terms of 

h k = i(u(x).n(x)) 

1 f dk dk ' 7 (u(k) • n(k')) 



2 



(2?r) d (2tt) c 
dkH K (k) 

Therefore, one dimensional magnetic helicity Hm is 

Ank 2 

HKik) = (^^(k). 
Using the definition fi(k) =ikx u(k), we obtain 



(m (k,t)tt i (k / ,«)) = 



P,j(k) C""(k) - - 



.2 



(27r) d <5(k + k'). 



Note that the magnetic helicity breaks mirror symmetry. 
We can Fourier transform in time as well using 

u (x, t) = J dku (k, ui) exp (ik • x — iuit) 
u(k, ui) — J dxdtu (x, i) exp (— ik • x + iuit) 
where dk — dkdw / (2n) d+1 . The resulting fluid equations in k = (k, ui) space are 



(-iu + vk 2 ) u t (k) = -^+ m (k) / dp [ Uj (p)u m (q)] , (25) 

After we have introduced the energy spectra and other second-order correlation functions, we move on to discuss 
their evolution. 



E. Energy Equations 

The energy equation for general (compressible) Navier-Stokes is quite complex. However, incompressible Navier- 
Stokes equation is relatively simpler, and is discussed below. 

From the evolution equations of fields, we can derive the following spectral evolution equations for incompressible 
MHD 

(!+**«) C(M - {d _ 1} | (k + v) / k , +p+q=o (Jp[-3 ((*' ■ u(q)) C(p) . uflO)) (26) 

where 9 stands for the imaginary part. Note that k' + p + q = and k' = — k. In Eq. J 2 (ill the first term in the RHS 
provides the energy transfer from the velocity modes to u(k) mode. Note that pressure couples with compressible 
modes, hence it is absent in the incompressible equations. 

In a finite box, using ^|u(k)| 2 ^) = C(k)/((rf — l)L d ), and <5(k)(27r) <i = L d (see Appendix A), we can show that 

(§i + 2vk2 ) \ (i u ( k )' 2 ) = E « k ' • u ^)) ( u (p) ■ u ( k ')))]' 

^ ' k'+p+q=0 

Many important quantities, e.g. energy fluxes, can be derived from the energy equations. We will discuss these 
quantities in the next section. 
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III. MODE-TO-MODE ENERGY TRANSFERS AND FLUXES IN MHD TURBULENCE 



In turbulence energy exchange takes place between various Fourier modes because of nonlinear interactions. Basic 
interactions in turbulence involves a wavenumber triad (k', p, q) satisfying k' + p + q = 0. Usually, energy gained by 
a mode in the triad is computed using the combined energy transfer from the other two modes [29]. Recently Dar et 
al. Q devised a new scheme to compute the energy transfer rate between two modes in a triad, and called it "mode- 
to-mode energy transfer". They computed energy cascade rates and energy transfer rates between two wavenumber 
shells using this scheme. We will review these ideas in this section. Note that we are considering only the interactions 
of incompressible modes. 



A. "Mode-to-Mode" Energy Transfer in Fluid Turbulence 



In this subsection we discuss evolution of energy in turbulent fluid in a periodic box. We consider ideal case where 
viscous dissipation is zero [y = 0). The equations are given in Lesieur j2^| 

~W (k')l 2 - E k > +P+q =o ~\* I(k' • u(q)) (u(k') • u(p)) + (k' • u(p)) (u(k') • u(q))] , (27) 

where 3 denotes the imaginary part. Note that the pressure does not appear in the energy equation because of the 
incompressibility condition. 

Consider a case in which only three modes u(k'), u(p), u(q), and their conjugates are nonzero. Then the above 
equation yields 

|l| U (k')| 2 = is««(k'|p,q), (28) 

where 

S(k'|p, q) = -3 [(k' • u(q)) (u(k') • u(p)) + (k' • u(p)) (u(k') • u(q))] . (29) 

Lesieur and other researchers physically interpreted S'(k'|p, q) as the combined energy transfer rate from modes p 
and q to mode k'. The evolution equations for |u(p)| 2 and |u(q)| 2 are similar to that for |w(k')| 2 . By adding the 
energy equations for all three modes, we obtain 



^[| U (k')| 2 + | U (p)| 2 + | U (q)| 2 ]/2 = 5 Ml (k'| P ,q) + ^"( P |q,k') + 5™(q|k',p) (30) 

= 3[(q.u(q))(u(k')-u(p)) (31) 
+ (p.u(p))(u(k')-u(q)) (32) 
+ (k'.u(k'))(u(p).u(q))] (33) 

For incompressible fluid the right-hand-side is identically zero because k' • u(k') = 0. Hence the energy in each 
interacting triad is conserved , i.e., 

|u (k')| 2 + \u (p)| 2 + \u (q)| 2 = const. (34) 

The question is whether we can derive an expression for mode-to-mode energy transfer rates from mode p to mode 
k', and from mode q to mode k' separately. Dar et al. @] showed that it is meaningful to talk about energy transfer 
rate between two modes. They derived an expression for the mode-to-mode energy transfer, and showed it to be 
unique apart from an irrelevant arbitrary constant. They referred to this quantity as "mode-to-mode energy transfer". 
Even though they talk about mode-to-mode transfer, they are still within the framework of triad interaction, i.e., a 
triad is still the fundamental entity of interaction. 



1. Definition of Mode-to-Mode Transfer in a Triad 



Consider a triad (k',p,q). Let the quantity i? Mtl (k'|p|q) denote the energy transferred from mode p to mode k' 
with mode q playing the role of a mediator. We wish to obtain an expression for R. 
The i?'s should satisfy the following relationships : 
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1. The sum of energy transfer from mode p to mode k' (i? uu (k'|p|q)), and from mode q to mode k' (i?™(k'|p|q)) 
should be equal to the total energy transferred to mode k' from modes p and q, i.e., «S" m (k'|p, q) [see Eq. ll28|l ]. 
That is, 

iT"(k'|p|q) + iT"(k'|q|p) = S""(k'|p, q), (35) 
iT"(p|k'|q) + i?™(p|q|k') = S""(p|k', q), (36) 

iT"(q|k'|p) + J R""(q|p|k') = 5™(q|k', P ). (37) 

2. By definition, the energy transferred from mode p to mode k', i? utl (k'|p|q), will be equal and opposite to the 
energy transferred from mode k' to mode p, i? mi (p|k'|q). Thus, 

iT u (k'|p|q) + iT u (p|k'|q) - 0, (38) 
iT u (k'|q|p) + i? ,m (q|k'|p) = 0, (39) 

i?" tl (p|q|k') + i? im (q| P |k') = 0. (40) 

These are six equations with six unknowns. However, the value of the determinant formed from the Eqs. (I35I40J) is 
zero. Therefore we cannot find unique i?'s given just these equations. In the following discussion we will study the 
set of solutions of the above equations. 

2. Solutions of equations of mode-to-mode transfer 

Consider a function 

S""(k'|p|q) = -S([k' • u(q)] [u(k') • u(p)]) (41) 

Note that ^""(k'lplq) is altogether different function compared to S^k'lp, q). In the expression for ^""(k'lplq), the 
field variables with the first and second arguments are dotted together, while the field variable with the third argument 
is dotted with the first argument. 

It is very easy to check that 5 uu (k'|p|q) satisfy the Eqs. 1135140)) . Note that 5""(k'|p|q) satisfy the Eqs. I|38l40j) 
because of incompressibility condition. The above results implies that the set of 5 u "(||)'s is one instance of the 
i? uli (||)'s. However, 5" u (k'|p|q) is not a unique solution. If another solution i?™(k'|p|q) differs from S'(k / |p|q) by an 
arbitrary function Xa, i-e., i?™(k'|p|q) = 5 UIl (k'|p|q) + Xa, then by inspection we can easily see that the solution 
of Eqs. i|35l40Jl must be of the form 

E""(k'|p|q) = 5 uu (k'|p|q)+Jf A (42) 
H utl (k'|q|p) = S w (k'|q|p) - X A (43) 
J R™(p|k'|q) = 5""(p|k'|q)-X A (44) 
i?™(p|q|k') = 5 wl (p|q|k')+^A (45) 
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R uu (k'lplq)=S uu (k'lplq)+X A 



Figure 1: Mode-to-mode energy transfer in fluid turbulence. <S™ u (k'|p|q) represents energy transfer rate from mode u(p) to 
mode u(k') with the mediation of mode u(q). Xa is the arbitrary circulating transfer. 



J? uu (q|k'|p) = S Mtl (q|k'|p) + X A (46) 



E""(q|p|k') = S uu (q|p|k')-X A (47) 

Hence, the solution differs from 5" u (k'|p|q) by a constant. See Fig. Q]for illustration. A careful observation of the 
figure indicates that the quantity Xa flows along p — > k' — > q — > p, circulating around the entire triad without 
changing the energy of any of the modes. Therefore we will call it the Circulating transfer. Of the total energy 
transfer between two modes, S uu + Xa, only S uu can bring about a change in modal energy. Xa transferred from 
mode p to mode k' is transferred back to mode p via mode q. Thus the energy that is effectively transferred from 
mode p to mode k' is just 5 u "(k'|p|q). Therefore 5 uu (k'|p|q) can be termed as the effective mode-to-mode energy 
transfer from mode p to mode k'. 

Note that Xa can be a function of wavenumbers k', p, q, and the Fourier components u(k'), u(p), u(q). It may be 
possible to determine Xa using constraints based on invariance, symmetries, etc. Dar et al. [a] attempted to obtain 
Xa using this approach, but could show that Xa is zero to linear order in the expansion. However, a general solution 
for Xa could not be found. Unfortunately, Xa cannot be calculated even by simulation or experiment, because we 
can experimentally compute only the energy transfer rate to a mode, which is a sum of two mode-to-mode energy 
transfers. The mode-to-mode energy transfer rate is really an abstract quantity, somewhat similar to "gauges" in 
electrodynamics. 

The terms UjdjUi and UiUjdjUi are nonlinear terms in the Navier-Stokes equation and the energy equation respec- 
tively. When we look at the formula (14 XB carefully, we find that the ttj(q) term is contracted with kj in the formula. 
Hence, Uj field is the mediator in the energy exchange between first (ttj) and third field (tt,) of UiUjdjUi. 

In this following discussion we will compute the energy fluxes and the shell-to-shell energy transfer rates using 
^"(k'lplq). 

B. Shell-to-Shell Energy Transfer in Fluid Turbulence Using Mode-to-mode Formalism 

In turbulence energy transfer takes place from one region of the wavenumber space to another region. Domaradzki 
and Rogallo have discussed the energy transfer between two shells using the combined energy transfer 5 uu (k'|p,q). 



12 




Figure 2: Shell-to-shell energy transfer from wavenumber-shell m to wavenumber-shell n. The triads involved in this process 
fall in two categories: Type I, where both p and q are inside shell m, and Type II, where only p is inside shell m. 



They interpret the quantity 

Tnm = ^EE S" u (k'| P) q). (48) 

k'6n p£m 

as the rate of energy transfer from shell to to shell n . Note that k'-sum is over shell n, p-sum over shell to, and 
k' + p + q = 0. However, Domaradzki and Rogallo ^3] themselves points out that it may not be entirely correct to 
interpret the formula lH8l) as the shell-to-shell energy transfer. The reason for this is as follows. 

In the energy transfer between two shells m and n, two types of wavenumber triads are involved, as shown in Fig. 
|2j The real energy transfer from the shell to to the shell n takes place through both k'-p and k'-q legs of triad I, but 
only through k'-p leg of triad II. But in Eq. l(4*8jl summation erroneously includes k'-q leg of triad II also along with 
the three legs given above. Hence Domaradzki and Ragallo's formalism pjj] do not yield totally correct shell-to-shell 
energy transfers, as was pointed out by Domaradzki and Rogallo themselves. We will show below how Dar et al.'s 
formalism Q overcomes this difficulty. 

By definition of the the mode-to-mode transfer function i? tt "(k'|p|q), the energy transfer from shell to to shell n 
can be defined as 

= EE ^ u ( k 'lp|q) (49) 

k'Gn p€m 

where the k'-sum is over the shell n, and p-sum is over the shell to. The quantity R uu can be written as a sum of an 
effective transfer ^""(k'lplq) and a circulating transfer X\. As discussed in the last section, the circulating transfer 
Xa does not contribute to the energy change of modes. From Figs. Hand El we can see that Xa flows from the shell 
to to the shell n and then flows back to to indirectly through the mode q. Therefore the effective energy transfer 
from the shell m to the shell n is just ^""(k'lplq) summed over all the k'-modes in the shell n and all the p-modes 
in the shell to, i.e., 

C=EE ^"(k'| P |q). (50) 

k'Gn p£m 
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Clearly, the energy transfer through k' — q of the triad II of Fig. [21 is not present in T"™ in Dar et al.'s formalism 
because q ^ to. Hence, the formalism of the mode-to-mode energy transfer rates provides us a correct and convenient 
method to compute the shell-to-shell energy transfer rates in fluid turbulence. 



C. Energy Cascade Rates in Fluid Turbulence Using Mode-to-mode Formalism 

The kinetic energy cascade rate (or flux) n in fluid turbulence is defined as the rate of loss of kinetic energy by 
the modes inside a sphere to the modes outside the sphere. Let k be the radius of the sphere under consideration. 
Kraichnan 22j, Leslie [^fJl, an d others have computed the energy flux in fluid turbulence using 5""(k'|p, q) 

E(*d) = - E E ^™(k'|p,q). (51) 

|k|<fc |p|>fco 

Although the energy cascade rate in fluid turbulence can be found by the above formula, the mode-to-mode approach 
of Dar et al. provides a more natural way of looking at the energy flux. Since i?""(k'|p|q) represents energy 
transfer from p to k' with q as a mediator, we may alternatively write the energy flux as 

n ( fc o) = E E ^(k'lplq). (52) 

|k'|>fc |p|<fe 

However, i? MU (k'|p|q) = ^""(k'lplq) + and the circulating transfer makes no contribution to the energy flux 
from the sphere because the energy lost from the sphere through returns to the sphere. Hence, 

n(*o) = E E ^(k'lplq)- (53) 

|k'|>feo |p|<fco 

Both the formulas given above, Eqs. l(5T|l and l(53"jl . are equivalent as shown by Dar et al. @. 
Frisch ^jj has derived a formula for energy flux as 

n(fco) = <u< • (u< • Vu> )) + <u< • (u> • Vu> )> . (54) 

It is easy to see that the above formula is consistent with mode-to-mode formalism. As discussed in the Subsection 
IIII A 21 the second field of both the terms are mediators in the energy transfer. Hence in mode-to-mode formalism, 
the above formula will translate to 

n (M = E E - s t( k ' ' u< (q)) ( u< (p) ' u> ( k ')) + ( k ' ' u> (q)) ( u< (p) ' u> ( k '))] < 

fe>feo p<fco 

which is same as mode-to-mode formula II53H of Dar et al. Q. 
The above quantities are computed numerically or theoretically. 



D. Digression to Infinite Box 

In the above discussion we assumed that the fluid is contained in a finite volume. In simulations, box size is typically 
taken to 2ir. However, most analytic calculations assume infinite box. It is quite easy to transform the equations 
given above to those for infinite box using the method described in Appendix. Here, the evolution of energy spectrum 
is given by (see Section Hill 

l + MJ ) C(k ' () " (^WTT?)/ k , +p+ ,. (^ [S ""< k ' |p| '' )l (55) 

The shell-to-shell energy transfer rate T nm from the TO-th shell to the n-th shell is 

r - = W+P+.,) L L <^ <s "" (k ' |p|q)> • (56) 

In terms of Fourier transform, the energy cascade rate from a sphere of radius fcois 

" wW+P + q) L <^F L <S "" <k ' IPlq)> ' <57) 
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For isotropic flows, after some manipulation and using Eq. (1241 . we obtain j2^| 

(^- t +2vkAE(k,t)=T{k,t), (58) 

where T(k,t), called transfer function, can be written in terms of S' YX (k'|p|q). The above formulas will be used in 
analytic calculations. 

The mode-to-mode formalism discussed here is quite general, and it can be applied to scalar turbulence p^ |. 
MHD turbulence, Rayleigh-Benard convection, enstrophy, Electron MHD etc. Some of these issues are discussed in 
Appendices C and D. One key assumption however is incompressibility. In the next section we will discuss various 
turbulence phenomenologies and models of fluid turbulence. 

IV. TURBULENCE PHENOMENOLOGICAL MODELS 

In the last two sections we introduced Navier-Stokes equation, and spectral quantities like the energy spectra and 
fluxes. These quantities have been analyzed using (a) phenomenological (b) numerical (c) analytical (d) experimental 
methods. In the present section we will present the most important phenomenological model called Kolmogorov's 
phenomenology of turbulence. 

A. Kolmogorov's 1941 Theory for Fluid Turbulence 

For ho mog eneous, isotropic, incompressible, and steady fluid turbulence with vanishing viscosity (large Re), Kol- 
mogorov [ijj, l2fl [2l|, [23] derived an exact relation that 

((Au)j) = ~el (59) 

where (Au)ii is component of u(x + 1) — u(x) along 1, e is the dissipation rate, and I lies between forcing scale (L) 
and dissipative scales (Id), i.e., « I « I. This intermediate range of scales is called inertial range. Note that the 
above relationship is universal, which holds independent of forcing and dissipative mechanisms, properties of fluid 
(viscosity), and initial conditions. Therefore it finds applications in wide spectrum of phenomena, e. g., atmosphere, 
ocean, channels, pipes, and astrophysical objects like stars, accretion disks etc. 

More popular than Eq. Q59I1 is its equivalent statement on energy spectrum. If we assume Au to be fractal, and e 
to be independent of scale, then 

((Auf) oc e 2 / 3 l 2 / 3 (60) 

Fourier transform of the above equation yields 

E(k) = K Ko e 2/3 k- r ° /3 (61) 

where Kko is a universal constant, commonly known as Kolmogorov's constant. Eq. l(6T|l has been supported by 
numerous experiments and numerical simulations. Kolmogorov's constant Kko has been found to lie between 1.4-1.6 
or so. It is quite amazing that complex interactions among fluid eddies in various different situations can be quite 
well approximated by Eq. 116 111 . 

Kolmogorov's derivation of Eq. ll59l) is quite involved. However, Eqs. l(59ll6T|i can be derived using scaling arguments 
(dimensional analysis) under the assumption that 

1. The energy spectrum in the inertial range does not depend on the large-scaling forcing processes and the small- 
scale dissipative processes, hence it must be a power law in the local wavenumber. 

2. The energy transfer in fluid turbulence is local in the wavenumber space. The energy supplied to the fluid at 
the forcing scale cascades to smaller scales, and so on. Under steady-state the energy cascade rate is constant 
in the wavenumber space, i. e., n(fc) = constant = e. 

In the framework of Kolmogorov's theory, several interesting deductions can be made. 
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1. Kolmogorov's theory assumes homogeneity and isotropy. In real flows, large-scales (forcing) as well as dissipative 
scales do not satisfy these properties. However, experiments and numerical simulations show that in the inertial 
range I), the fluid flows are typically homogeneous and isotropic. 

2. The velocity fluctuations at any scale I goes as 

ui ra e 1 / 3 / 1 / 3 . (62) 
Therefore, the effective time-scale for the interaction among eddies of size I is 

Tl » - « e~ l W 3 . (63) 
ui 



3. An extrapolation of Kolmogorov's scaling to the forcing and the dissipative scales yields 



Taking f w u; d /d, one gets 

,3\ V4 



c 



*d « ( - I • (65) 



Note that the dissipation scale, also known as Kolmogorov's scale, depends on the large-scale quantity e apart 
from kinematic viscosity. 

4. From the definition of Reynolds number 

U L L U L L fL\ 4/3 . 
Re = « « ( - ) (66) 

Therefore, 



r « i?e 3 / 4 . (67) 
<d 

Onset of turbulence depends on geometry, initial conditions, noise etc. Still, in most experiments turbulences 
sets in after Re of 2000 or more. Therefore, in three dimensions, number of active modes (L/ld) 3 is larger than 
26 million. These large number of modes make the problem quite complex and intractable. 

Space dimension does not appear in the scaling arguments. Hence, one may expect Kolmogorov's scaling to 
hold in all dimensions. It is however found that the above scaling law is applicable in three dimension only. In 
two dimension (2D), conservation of enstrophy changes the behaviour significantly (see next two sections). The 
solution for one-dimensional incompressible Navier-Stokes is u(x, t) = const, which is a trivial solution. 

Mode-to-mode energy transfer term S(k\p\q) measures the strength of nonlinear interaction. Kolmogorov's 
theory implicitly assumes that energy cascades from larger to smaller scales. It is called local energy transfer in 
Fourier space. These issues will be discussed in Section fVIIII 

Careful experiments show that the spectral index is close to 1.71 instead of 1.67. This correction of « 0.04 
is universal and is due to the small-scale structures. This phenomena is known as intermittency, and will be 
discussed in Section HXl 



8. Kolmogorov's model for turbulence works only for incompressible flow. It is connected to the fact that incom- 
pressible flow has local energy transfer in wavenumber space. Note that Burgers equation, which represents 
compressible flow (U C s ), has k~ 2 energy spectrum, very different from Kolmogorov's spectrum. 

Kolmogorov's theory of turbulence had a major impact on turbulence research because of its universality. Properties 
of scalar, MHD, Burgers, Electron MHD, wave turbulence have been studied using similar arguments. 

As discussed in earlier sections, apart from energy spectra, there are many other quantities of interest in turbulence. 
Some of them are kinetic helicity, enstrophy etc. The statistical properties of these quantities are quite interesting, 
and they are addressed using Absolute Equilibrium State discussed below. 
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B. Absolute Equilibrium States 

In fluid turbulence when viscosity is identically zero (inviscid limit), kinetic energy is conserved in the incompressible 
limit. Now consider independent Fourier modes (transverse to wavenumbers) as state variables y a (t). Lesieur [2<| has 
shown that these variables move in a constant energy surface, and the motion is area preserving like in Liouville's 
theorem. Now we look for equilibrium probability-distribution function P{{y a }) for these state variables. Once we 
assume ergodicity, the ideal incompressible fluid turbulence can be mapped to equilibrium statistical mechanics [2fl| . 

By applying the usual arguments of equilibrium statistical mechanics we can deduce that at equilibrium, the 
probability distribution function will be 

P(y 1 ,...,y m ) = iexp ^-^ct^^Y 

where a is a positive constant. The parameter a corresponds to inverse temperature in the Boltzmann distribution. 
Clearly 

(yl) = J IkdyiylPiivi}) = 1, 

independent of a. Hence energy spectrum C(k) is constant, and 1-d spectrum will be proportional to /c d_1 j2j|. This 
is very different from Kolmogorov's spectrum for large Re turbulence. Hence, the physics of turbulence at v = 
(inviscid) differs greatly from the physics at v — > 0. This is not surprising because (a) turbulence is a nonequilibrium 
process, and (b) Navier-Stokes equation is singular in v. 

Even though nature of inviscid flow is very different from turbulent flow, Kraichnan and Chen [26] suggested that 
the tendency of the energy cascade in turbulent flow could be anticipated from the absolute equilibrium states. Using 
absolute equilibrium theory, Kraichnan j^l showed that in two dimensions, enstrophy cascades forward, but energy 
cascades backward (see also Lesieur ^3)- The above prediction holds good for real fluids. 



V. EXPERIMENTAL RESULTS ON TURBULENCE 



Analytical results are very rare in turbulence research because of complex nature of turbulence. Therefore, experi- 
ments and numerical simulations play very important role in turbulence research. In fluid turbulence, engineers have 
been able to obtain necessary information from experiments (e.g., wind tunnels), and successfully design complex 
machines like aeroplanes, spacecrafts etc. This aspect of fluid turbulence is not being covered here. For details on 
experiments, refer to books on fluid turbulence, e.g., Davidson @|. 
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VI. NUMERICAL INVESTIGATION OF FLUID TURBULENCE 



Like experiments, numerical simulations help us test existing models and theories, and inspire new one. In addition, 
numerical simulations can be performed for conditions which may be impossible in real experiments, and all the field 
components can be probed everywhere, and at all times. Recent exponential growth in computing power has fueled 
major growth in this area of research. Of course, numerical simulations have limitations as well. Even the best 
computers of today cannot resolve all the scales in a turbulent flow. We will investigate these issues in this section. 

There are many numerical methods to simulate turbulence on a computer. Engineers have devised many clever 
schemes to simulate flows in complex geometries; however, their attention is typically at large scales. Physicists 
normally focus on intermediate and small scales in a simple geometry because these scales obey universal laws. Since 
nonlinear equations are generally quite sensitive, one needs to compute both the spatial and temporal derivatives 
as accurately as possible. It has been shown that spatial derivative could be computed "exactly" using Fourier 
transforms given enough resolutions Q. Therefore, physicists typically choose spectral method to simulate turbulence. 
Note however that several researchers have used higher order finite-difference scheme and have obtained comparable 
results. 



A. Numerical Solution of fluid Equations using Pseudo-Spectral Method 

In this subsection we will briefly sketch the spectral method for 3D flows. For details refer to Canuto et al. Q|. 
The fluid equations in Fourier space is written as 

d U ip^ = -zkp(k,t)-FT[u(k,t)-Vu(k,t)]-^fc 2 u(k,t) + f 

where FT stands for Fourier transform, and f(k, t) is the forcing function. The flow is assumed to be incompressible, 
i. e., k- u(k, t) — 0. We assume periodic boundary condition with real-space box size as (2ir) x (2ir) x (2ir), and 
Fourier-space box size as (nx, ny, nz). The allowed wavenumbers are k = (k x , k y , k z ) with k x — (~n x /2 : n x /2), k y = 
( — n,j/2 : it „/•>). /.-.. = ( — a-; /2 : n «/^- The reality condition implies that z ( -k) z ! » (k), therefore, we need 



to consider only half of the modes Typically we take {—n x /2 : n x /2, —n y /2 : n y /2,0 : n z /2), hence, we have 
N — n x * n y * (n z /2 + 1) coupled ordinary differential equations. The objective is to solve for the field variables at a 
later time given initial conditions. The following important issues are involved in this method: 

1. The Navier-Stokes equation is converted to nondimensionalized form, and then solved numerically. The param- 
eter v is inverse Reynold's number. Hence, for turbulent flows, v is chosen to be quite small (typically 1CP 3 or 
1CP 4 ). In Section liV Al we deduced using Kolmogorov's phenomenology that the number of active modes are 

If we choose a moderate Reynolds number v~ x = 10 4 , N will be 10 9 , which is a very large number even for 
the most powerful supercomputers. To overcome this difficulty, researchers apply some tricks; the most popular 
among them are introduction of hyperviscosity and hyperresistivity, and large-eddy simulations. Hyperviscous 
(hyperresistive) terms are of the form (vj)k J u (k) with j > 2; these terms become active only at large wavenum- 
bers, and are expected not to affect the inertial range physics, which is of interest to us. Because of this property, 
the usage of hyperviscosity and hyperresistivity has become very popular in turbulence simulations. Large-eddy 
simulations are discussed in various books (e.g., see Pope [4l| . Just to note, one of the highest resolution fluid 
turbulence simulation is by Gotoh [l£( on a 4096 3 grid; this simulation was done on Fujitsu VPP5000/56 with 
32 processors with 8 Gigabytes of RAM on each processor, and it took 500 hours of computer time. 

2. The computation of the nonlinear terms is the most expensive part of turbulence simulation. A naive calculation 
involving convolution will take 0(N 2 ) floating point operations. It is instead efficiently computed using Fast 
Fourier Transform (FFT) as follows: 

(a) Compute u (x) from u (k) using Inverse FFT. 

(b) Compute Ui(x)uj(x) in real space by multiplying the fields at each space points. 

(c) Compute FFT[ui(x)u.,-(x)] using FFT. 

(d) Compute ikjFFT[ui(~x)uj{'x)] by multiplying by kj and summing over all j. This vector is 
-FFT [u (k, t) ■ Vu (k, t)] . 

Since FFT takes 0(N log N), the above method is quite efficient. The multiplication is done in real space, 
therefore this method is called pseudo-spectral method instead of just spectral method. 
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3. Products itj(x)itj(x) produce modes with wavenumbers larger than k max . On FFT, these modes get aliased 
with k < k max and will provide incorrect value for the convolution. To overcome this difficulty, last 1/3 modes 
of fields z i (k) are set to zero (zero padding), and then FFTs are performed. This scheme is called 2/3 rule. For 
details refer to Canuto et al. 

4. Pressure is computed by taking the dot product of Navier-Stokes equation with k. Using incompressibility 
condition one obtains 

p(k,t) = ^-FT[u(x I t).Vu(x I t)]. 

To compute p(k) we use already computed nonlinear term. 

5. Once the right-hand side of the Navier-Stokes equation could be computed, we could time advance the equation 
using one of the standard techniques. The viscous terms are advanced using an implicit method called Crank- 
Nicholson's scheme. However, the nonlinear terms are advanced using Adam-Bashforth or Runge-Kutta scheme. 
One uses either second or third order scheme. Choice of dt is determined by CFL criteria (dt < (Ax)/U rms ). 
By repeated application of time-advancing, we can reach the desired final time. 

6. When forcing f = 0, the total energy gets dissipated due to viscosity. This is called decaying simulation. On 
the contrary, forced simulation have nonzero forcing (f ^ 0), which feed energy into the system, and the system 
typically reaches a steady-state in several eddy turnover time. Forcing in turbulent systems are typically at large- 
scale eddies (shaking, stirring etc.). Therefore, in forced turbulence f is typically applied at small wavenumbers, 
which could feed kinetic energy and kinetic helicity. 

Spectral method has several disadvantages as well. This method can not be easily applied to nonperiodic flows. That 
is the reason why engineers hardly use spectral method. Note however that even in aperiodic flows with complex 
boundaries, the flows at small length-scale can be quite homogeneous, and can be simulated using spectral method. 
Spectral simulations are very popular among physicists who try to probe universal features of small-scale turbulent 
flows. 

The numerical results on turbulent energy spectrum and fluxes are described in many turbulence literature, e. g., 
Lesieur j2^|. In Section 8 we describe a numerical result on shell-to-shell energy transfer in 512 3 simulation. In recent 
times a technique called large-eddy simulation (LES) has become very popular. LES enables us to perform turbulence 
simulations on smaller grids. In this paper we donot cover this topic. 

In the next three sections we will describe the field-theoretic calculation of renormalized viscosity. 
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VII. RENORMALIZATION GROUP ANALYSIS OF FLUID TURBULENCE 

In Section ?? we discussed various existing turbulence models. Here we will describe some field-theoretic calcula- 
tions. 

Field theory is well developed, and has been applied to many areas of physics, e.g., Quantum Electrodynamics, 
Condensed Matter Physics etc. In this theory, the equations are expanded perturbatively in terms of nonlinear term, 
which are considered small. In fluid turbulence the nonlinear term is not small; the ratio of nonlinear to linear 
(viscous) term is Reynolds numbers, which is large in turbulence regime. This problem appears in many areas of 
physics including Quantum Chromodynamics (QCD), Strongly Correlated Systems, Quantum Gravity etc., and is 
largely unsolved. To overcome the above difficulty, some clever schemes have been adopted such as Direct Interaction 
Approximation, Renormalization Groups (RG), Eddy-damped quasi-normal Markovian approximations, etc. We 
discuss some of them below. A simple-minded calculation of Green's function shows divergence at small wavenumbers 
(infrared divergence). One way to solve this problem is by introducing an infrared cutoff for the integral. The reader 
is referred to Leslie [j^l f° r details. RG technique, to be described below, is a systematic procedure to cure this 
problem. 

A. Renormalization Groups in Turbulence 

Renormalization Group Theory (RG) is a technique which is applied to complex problems involving many length 
scales. Many researchers have applied RG to turbulence. Over the years, several different RG applications for 
turbulence has been discovered. Broadly speaking, they fall in three different categories: 

Yakhot-Orszag (YO) Perturbative approach 

Yakhot and Orszag's [58] work, motivated by Forster et al. ^3] an d Fournier and Frisch ^4[> ls the first compre- 
hensive application of RG to turbulence. It is based on Wilson's shell-elimination procedure. Also refer to Smith and 
Woodruff [44] for details. Here the renormalized parameter is function of forcing noise spectrum D(k) — D k~ y . It 
is shown that the local Reynolds number A is 



^ 3 (A)A £ ' 

where Ao is the expansion parameter, A is the cutoff wavenumber, and e = 4 + y — d [58j. It is found that ^(A) increases 
as A decreases, therefore, A remains small (may not be less that one though) compared to Re as the wavenumber 
shells are eliminated. Hence, the "effective" expansion parameter is small even when the Reynolds number may be 
large. 

The RG analysis of Yakhot and Orszag [[jjjj yielded Kolmogorov's constant Kk — 1.617, turbulent Prandtl number 
for high-Reynolds-number heat transfer Pt = 0.7179, Batchelor constant Ba = 1.161 etc. These numbers are quite 
close to the experimental results. Hence, Yakhot and Orszag's method appears to be highly successful. However there 
are several criticisms to the YO scheme. Kolmogorov's spectrum results in the YO scheme for e = 4, far away from 
e = 0, hence epsilon-expansion is questionable. YO proposed that higher order nonlinearities are "irrelevant" in the 
RG sense for e = 0, and are marginal when e = 4. Eyink objected to this claim and demonstrated that the higher 
order nonlinearities are marginal regardless of e. Kraichnan (23 compared YO's procedure with Kraichnan's Direct 
Interaction Approximation [231 and raised certain objections regarding distant-interaction in YO scheme. For details 
refer to Zhou and McComb |f[jj and Smith and Woodruff 0. 

Self- consistent approach of McComb and Zhou 

This is one of the nonperturbative method, which is often used in Quantum Field theory. In this method, a self- 
consistent equation of the full propagator is written in terms of itself and the proper vertex part. The equation may 
contain many (possibly infinite) terms, but it is truncated at some order. Then the equation is solved iteratively. 
McComb 33], Zhou and coworkers j^| have applied this scheme to fluid turbulence, and have calculated renormalized 
viscosity and Kolmogorov's constant successfully. Direct Interaction Approximation of Kraichnan is quite similar to 
self-consistent theory (Smith and Woodruff 0]). 

The difficulty with this method is that it is not rigorous. In McComb and Zhou's procedures, the vertex correction 
is not taken into account. Verma (49I IHnL l5ll| has applied the self-consistent theory to MHD turbulence. 
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Figure 3: The wavenumber shells to be averaged during renormalization procedure. 



Callan-Symanzik Equation for Turbulence 

DeDominicis and Martin @ and Teodorovich ^3] obtained the RG equation using functional integral. Teodorovich 
obtained Kko = 2.447, which is in not in good agreement with the experimental data, though it is not too far away. 
It has been shown that Wilson's shell-renormalization and RG through Callan-Symanzik equation are equivalent 
procedure. However, careful comparison of RG schemes in turbulence is not completely worked out. 

In the following discussion we will discuss McComb's RG scheme is some detail. The other schemes have been 
discussed in great lengths in several books and review articles. After renormalization, in Section IVIIII we will discuss 
the computation of energy fluxes. These calculations are done using self-consistent field theory, a scheme very similar 
to DIA. At the end we will describe Eddy-damped quasi-normal Markovian approximation, which is very similar to 
the energy flux calculation. 

B. Physical Meaning of Renormalization in Turbulence 

The field theorists have been using renormalization techniques since 1940s. However, the physical meaning of renor- 
malization became clear after path-breaking work of Wilson [^3] ■ Here renormalization is a variation of parameters 
as we go from one length scale to the next. Following Wilson, renormalized viscosity and resistivity can also be 
interpreted as scale-dependent parameters. We coarse-grain the physical space and look for an effective theory at 
a larger scale. In this method, we sum up all the interactions at smaller scales, and as a outcome we obtain terms 
that can be treated as a correction to viscosity and resistivity. The corrected viscosity and resistivity are called 
"effective" or renormalized dissipative parameters. This procedure of coarse graining is also called shell elimination 
in wavenumber space. We carry on with this averaging process till we reach inertial range. In the inertial range the 
"effective" or renormalized parameters follow a universal powerlaw, e. g., renormalized viscosity u(l) oc I 4 / 3 . This is 
the renormalization procedure in turbulence. Note that the renormalized parameters are independent of microscopic 
viscosity or resistivity. 

In viscosity renormalization the large wavenumber shells are eliminated, and the interaction involving these shells 
are summed. Hence, we move from larger wavenumbers to smaller wavenumbers. However, it is also possible to go 
from smaller wavenumbers to larger wavenumber by summing the smaller wavenumber shells, e.g, for shear flows. 
This process is not coarse-graining, but it is a perfectly valid RG procedure, and is useful when the small wavenumber 
modes (large length scales) are linear. This scheme is followed in Quantum Electrodynamics (QED), where the 
electromagnetic field is negligible at a large distance (small wavenumbers) from a charge particle, while the field 
becomes nonzero at short distances (large wavenumber). In QED, the charge of a particle gets renormalized when 
we come closer to the charge particle, i. e., from smaller wavenumbers to larger wavenumbers. See Fig. OH for an 
illustration of wavenumber shells to be averaged. In the following subsection we will calculate renormalized viscosity 
using RG procedure. 
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C. Renormalization of viscosity using self-consistent procedure 

In this subsection we compute renormalized viscosity using self-consistent procedure. This work was done by 
McComb and his group workers. The renormalization of viscosity is performed from large wavenumber to smaller 
wavenumbers. 

McComb and his group workers took the following form of Kolmogorov's spectrum for kinetic energy 

E(k) = K Ko n 2 ^k- 5 / 3 , (69) 

where Kko is Kolmogorov's constant, and II is the total energy flux. The incompressible fluid equations in the Fourier 
space are 

(-iu + vk 2 )uAk) = -Ji# m (k) / dp[ Uj {p)u m {q)), (70) 

where 

P+ m (k) = kjPi^ + kmPi^k), (71) 

Here v is viscosity, and d is the space dimensionality. 

In this RG procedure the wavenumber range (k]y,ko) is divided logarithmically into N shells. The nth shell is 
(k n ,k„-i) where k n = h n ko (h < 1). In the following discussion, the elimination of the first shell (fei,fco) is carried 
out, and modified NS equation is obtained. Then one proceeds iteratively to eliminate higher shells and get a general 
expression for the modified fluid equation. The renormalization group procedure is as follows: 

1. The the spectral space is divided in two parts: 1. the shell (ki,ko) — k> , which is to be eliminated; 2. 
(kN,ki) — /c < , set of modes to be retained. Note that V(q) denote the viscosity and resistivity before the 
elimination of the first shell. 

2. Rewrite Eqs. lT7fil) for /c < and fc > . The equations for uf(k) and bf (k) modes are 

(-iw + S (0) (*))<(*) = - l -P+ m (k) J dp([u<(p)u<(k~p)} 

+2[uf[p)u>{k-p)] + {u>(p)u>(k~p)}) (72) 

The Ss appearing in the equations are usually called the "self-energy" in Quantum field theory language. In the 
first iteration, S( ) = V(o)k 2 . The equation for uf(k) modes can be obtained by interchanging < and > in the 
above equations. 

3. The terms given in the second and third brackets in the Right-hand side of Eqs. I|72J1 are calculated perturbatively. 
Since we are interested in the statistical properties of u fluctuations, we perform the usual ensemble average of 
the system j^. It is assumed that u > (fc) has gaussian distributions with zero mean, while u(fe) is unaffected 
by the averaging process. Hence, 

(u>{k)) = (73) 
(uf{k)) = «<(*) (74) 

and 

(u>(p)u>(q)) - P l3 (p)C(p)5(p + q) (75) 

The triple order correlations (uf (k)u> (p)u^ l (q)j are zero due to Gaussian nature of the fluctuations. Here, X 

stands for u or b. In addition, we also neglect the contribution from the triple nonlinearity ^u < (fc)u^(p)u^(g)^, 

as done in many of the turbulence RG calculations (33l IHsj. The effects of triple nonlinearity can be included 
following the scheme of Zhou and Vahala [fi5l |. 
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4. To the first order, the second bracketed terms of Eqs. I|72JI vanish, but the nonvanishing third bracketed terms 
yield corrections to Es. Refer to Appendix C for details. Eqs. lT72)l can now be approximated by 

(-w + S (0) +<SS (0) )u,<(fe) - -^+„(k) J dp[uf(p)u<(k-p)] (76) 

with 

1 ' A 



SE$(k) = -— — / dp[S(k,p,q)G(p)Cm (77) 
where 

S(k,p, q) = kp((d - 3)z + 2z 3 + (d- l)xy). (78) 
The integral A is to be done over the first shell. 

5. The frequency dependence of the correlation function is taken as: C(k,ui) — 2C(k)dt(G(k,uj)). In other words, 
the relaxation time-scale of correlation function is assumed to be the same as that of corresponding Green's 
function. Since we are interested in the large time-scale behaviour of turbulence, we take the limit uj going to 
zero. Under these assumptions, the frequency integration of the above equations yield 

r rM _ 1 f A dp S(k,p,q)C(q) 

° V ^ } (d - 1)P 7 p+q=k (2kY v [Q) {p)p* + V(0) {q)q* [n)) 

Note that v(k) = Y^ uu (k)/k 2 . There are some important points to remember in the above step. The frequency 
integral in the above is done using contour integral. It can be shown that the integrals are nonzero only when 
both the components appearing the denominator are of the same sign. For example, first term of Eq. lf79l) is 
nonzero only when both vto){p) an d J / (o)(9) are of the same sign. 

6. Let us denote f(i)(fc) as the renormalized viscosity after the first step of wavenumber elimination. Hence, 

f(i)(fc) = f(0)(*) + *f(0)(*); (80) 
We keep eliminating the shells one after the other by the above procedure. After n + 1 iterations we obtain 

" (n+i) ( fc ) = v {n) (k) + 8v {n) (k) (81) 

where the equation for 8v^ n ){k) is the same as the Eqs. lT79ll except that V(o){k) appearing in the equation is to 
be replaced by V( n )(k). Clearly i>( n+ i)(fc) is the renormalized viscosity and resistivity after the elimination of 
the (n + l)th shell. 

7. We need to compute <5f( n ) for various n. These computations, however, require vt n ). In our scheme we solve 
these equations iteratively. In Eqs. ifTfll ??) we substitute C{k) by one dimensional energy spectrum E{k) 

where Sd is the surface area of d-dimensional spheres. We assume that E{k) follows Eqs. . Regarding i/r n ), 
we attempt the following form of solution 

v {n) {k n k') = (^ ) 1/2 n 1 / 3 fc- 4 / 3 ^ n) (fc') 

with k — kn+ik' (k' < 1). We expect v* n Jk') to be a universal functions for large n. The substitution of 
C(fc), i/(„)(fc) yields the following equations: 

W) - (j~T) L,.^^ |S|tV ' fl . ;j (W^ k ^(W 1 (82) 

"(»+«(*') = hV^W + h-VH^ik') (83) 

where the integrals in the above equations are performed iteratively over a region 1 < p',q' < l/h with the 
constraint that p' + q' = k'. Fournier and Frisch 0| showed the above volume integral in d dimension to be 

J _^dp' = S d -iJdp'dq'(^pj (sina) d - 3 , (84) 
where a is the angle between vectors p' and q'. 
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Figure 4: Plot of u*(k') vs. k' for 2D and 3D fluid turbulence. In 2D, v is negative. 
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Figure 5: The plot of asymptotic ^*vs. d . For large d, the plot fits quite well with predicted d 1//2 curve. Adopted from 
Verma |H3|. 

8. Now the above equations are solved self-consistently with h = 0.7. This value is about middle of the range 
(0.55-0.75) estimated to be the reasonable values of h by Zhou et al. One starts with constant value of v* y 
and compute the integrals using Gauss quadrature technique. Once Su^ has been computed, v*^ is computed. 
This process is iterated till vt m+l Ak') ~ v* m ^(k'), that is, till they converge. The result of our RG analysis is 
given below. 

McComb and coworkers (IE HE EH successfully applied the above self-consistent renormalization group theory to 
2D and 3D fluid turbulence. They found that v*(k') converges quite quickly. For 3D the value of v*(k' — » 0) is 
approximately 0.38. See Fig. 01 for an illustration. 

For 2D turbulence v*(k') is negative as shown in Fig. Q] The function v* is not very well behaved as k' — > 0. Still, 
negative renormalized viscosity is consistent with negative eddy viscosity obtained using Test Field Model [24[ and 
EDQNM calculations We estimate v* w -0.60 . 

For large d , v* = rj*, and it decreases as dr 1 / 2 (see Fig. EJ; v* for pure fluid turbulence also decreases as g?" 1 / 2 , as 
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shown in the same figure. This is evident from Eqs. I|82H lsj. For large d 



d p>dq'(^Pj d * (sma) d - 3 ... ~ (T 1 / 2 , (85) 
S d -i 1 ( d^ 1 ' 2 



which leads to 



{d-l) 2 S d d?\2ir J 
S,-S 6 ,-Ss,Sg(k',p',q') = kpd(z + xy), (86) 



-, / j N 1/2 

or \2tt / 



hence is* ct d 1 I 2 . 

In conclusion the above RG procedure shows that 



£(fc) = K Ko Tl 2 / 3 k- 5 / 3 , (87) 

v{k = k n k') = K^n^k-^v'ik'), 



is a consistent solution of renormalization group equation. Here, Kko is Kolmogorov's constant, II is the energy flux, 
and v*(k') is a universal function that is a constant as k! — > 0. 



1. Helical turbulence 



Helical turbulence is defined for space dimension d = 3. We can extend the above the RG analysis to helical 
turbulence (Zhou |8l]). All the steps are the same except Eqs. lf75)) are replaced by 



(u>(p)u>(q)) = 



{2n)H(p + q) (89) 



Because of helicities, the equation for change in renormalized self-energy H79fl gets altered to 

c m = 1 f A dp S(k,p,q)C uu (q) + S'(k,p,q)H K (q) 

" (0)( > (d 1)P i p+q=k (2^ 1 v {0) (p)pi + m (q)? 

where defined below can be shown to be zero. 

S'(k,p,q) = P b + jm (k)P+ ab (p)e jal q l =0, 

The argument for vanishing of S' is follows. Since 8v is a proper scalar, and Hk is a pseudo scalar, S'(k,p,q) will 
be also be a pseudo scalar. In addition, S'(k,p,q) are also linear in k,p and q. This implies that Sl(k,p,q) must 
be proportional to q • (k x p), which will be zero because k = p + q. Hence S'(k,p, q) turn out to be zero. Hence, 
helicities do not alter the already calculated S(u)/ n \(k) in the earlier section. 

In fluid turbulence, there are some other interesting variations of field-theoretic calculations by DeDominicis and 
Martin jfl], Bhattacharjee 0, Carati @ and others. 

In the next section we will compute energy fluxes for fluid turbulence using field-theoretic techniques. 
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VIII. FIELD-THEORETIC CALCULATION OF ENERGY FLUXES AND SHELL-TO-SHELL ENERGY 

TRANSFER 

In this section we present calculation of energy flux using field-theoretic method. We assume the turbulence to be 
homogeneous and isotropic. Even though the real-world turbulence do not satisfy these properties, many conclusions 
drawn using these assumption provide us with important insights into the energy transfer mechanisms at small scales. 
The field-theoretic procedure requires Fourier space integrations of functions involving products of energy spectrum 
and the Greens functions. Since there is a general agreement on Kolmogorov-like spectrum for fluid turbulence, 
E(k) oc fc -5 / 3 is taken for the energy spectrum. For the Greens function, we substitute the "renormalized" or 
"dressed" Greens function computed in the previous section 

|U (see Section IVTTCl 



A. Calculation of Energy Flux 

As described in Section IIIII the energy flux from a wavenumber sphere of radius fco to the outside of the sphere of 
the same radius is 

n *> - (2.)W + p +q) L„ L is "" ik ' lvM) {m 

We assume that the kinetic energy is forced at small wavenumbers. 

We analytically calculate the above energy fluxes in the inertial range to leading order in perturbation series. It was 
assumed that u(k) is quasi-gaussian as in EDQNM approximation. Under this approximation, the triple correlation 
{XXX) is zero to zeroth order, but nonzero to first oder. To first order (XXX) is written in terms of (XX XX), 
which is replaced by its Gaussian value, a sum of products of second-order moment. Consequently, the ensemble 
average of S, (S), is zero to the zeroth order, but is nonzero to the first order. The first order terms for (S(k\p\q)) 
in terms of Feynman diagrams are given in Appendix C. They are given below in terms of Green's functions and 
correlation functions: 



(S(k\p\q)) = [ dt , (27r) d [T 1 (k,p,q)G(k,t-t')C(p,t,t')C(q,t,t') 

J —OO 

+T 2 (k, Pl q)G(p, t - t')C(k, t, t')C(q, t, t') 

+T 3 (k,p, q)G{q, t - t')C(k, t, t')C(p, t, t')]S(k' + p + q) (91) 

where Ti(k,p, q) are functions of wavevectors k,p, and q given in Appendix B. 

The Greens functions can be written in terms of "effective" or "renormalized" viscosity v(k) and resistivity r/(k) 
computed in Section IVlII 

G(k, t-t') = 0(t - t') exp (~v(k)k 2 (t - t')) 

The relaxation time for C(k,t,t') is assumed to be the same as that of G(k,t,t'). Therefore the time dependence of 
the unequal-time correlation functions will be 

C uu {k,t,t') = 6(t-t')eyLv(~v{k)k 2 (t-t')) C uuM (k,t,t) 

The above forms of Green's and correlation functions are substituted in the expression of (S), and the t' integral is 
performed. Now Eq. I|9()|) yields the following flux formula for IT(£;o): 

n( *» ) - L L w? - w* + + g x |r ' (t ' p - " c ""' c ""» 

+T 2 (k,p,q)C uu (k)C uu (q)+T 3 (k 1 p,q)C uu (k)C uu (p)}. (92) 

The expressions for the other fluxes can be obtained similarly. 

The equal-time correlation function C(k,t,t) at the steady-state can be written in terms of one dimensional energy 
spectrum as 

G^^^-^Eik), 
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where Sd is the surface area of d-dimensional unit spheres. We are interested in the fluxes in the inertial range. There- 
fore, we substitute Kolmogorov's spectrum [Ens. (l69?l ] for the energy spectrum. The effective viscosity is proportional 
to fc- 4 / 3 , i.e., 

v{k) = (K u ) x l 2 rtl z k^l z v\ (93) 

and the parameter v* was calculated in Section IVl!! 
We nondimensionalize Eq. 1(92)) by substituting jsFj] 

U ka k ° k ° /nA\ 

k = — : p = — v; q = — w. (94) 

u u u 



Application of Eq. lIMjl yields 



n = (K Ko f/ 2 u 

where the integral F(v,w) are 



(d~iys d j 



1 rl+v 

dv\n(l/v) / dw(vw) d ~ 2 (sin a) d ~ 3 F(v,w) 

Jl-V 



(95) 



F = ^TI 273 , ^ [ti(v,w)(vw)- d -% +h(v,w)w- d -i +t 3 (v,w)v- d -i], (96) 

Here ti(v, w) — Ti(k, kv, kw)/k 2 . Note that the energy fluxes are constant, consistent with the Kolmogorov's picture. 
We compute the bracketed term (denoted by I) numerically using Gaussian-quadrature method, and found it to be 
convergent. Using / the constant Kko can be calculated as 

K Ko = (iy 2/3 . (97) 

For 3D turbulence, the value of constant Kko computed using Eqs. I|95l I97|l is 1.58. This number is very good 
agreement with numerical and experimental estimate of Kolmogorov's constant. For 2D turbulence, we substitute 
v* = —0.60 in the above equations. The computation yields K\P w 6.3. 

For large d 

dp'dq'f^-Y 2 (sinaf- 3 ... ~ d' 1 ' 2 , (98) 



Sd-i J- 

d? V 2tt 



(d - l) 2 S d 



1/2 



v* = ii* - d- 1 ' 2 (99) 
t x = -h = kpd(z + xy), (100) 

and ^3 = 0. Using Eq. IIOOII and by matching the dimensions, it can be shown that K oc d -1 / 3 . This result is due to 
Fournier et al. [13] . 

All the above conclusions are for large Reynolds number or v — > limit. The behaviour of Navier-Stokes equation 
for viscosity v = (inviscid) is very different, and has been analyzed using absolute equilibrium theory (see Section 
IIV Bll . It can be shown using this theory that under steady state, energy is equipartitioned among all the modes, 
hence C(k) = const j^|. Using this result we can compute mode-to-mode energy transfer rates (S uu (k\p\q)) to first 
order in perturbation theory (Eq. [S]), which yields 

oc / CTi(fc,P, «) +Jf>P> g) + T 9 (k,P, q)) Const = 
J v(k)k 2 + v(p)p 2 + v{q)q 2 

because T\{k,p, q) +T 5 (k,p, q) +T 9 (k,p, q) = 0. Hence, under steady-state, their is no energy transfer among Fourier 
modes in inviscid Navier-Stokes. In other words "principle of detailed balance" holds here. Note that the above result 
holds for all space dimensions. Contrast this result with the turbulence situation when energy preferentially gets 
transferred from smaller wavenumber to larger wavenumber. This example contrasts equilibrium and nonequilibrium 
systems. 

After completing the discussion on energy fluxes for fluid turbulence, we now move on to theoretical computation 
of shell-to-shell energy transfer. 
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B. Field-theoretic Calculation of Shell-to-shell Energy Transfer 



Energy transfers between wavenumber shells provide us with important insights into the dynamics of turbulence. 
Kolmogorov's fluid turbulence model is based on local energy transfer between wavenumber shells. There are several 
quantitative theories in fluid turbulence about the amount of energy transfer between neighbouring wavenumber 
shells. For examples, Kraichnan |24| showed that 35% of the energy transfer comes from wavenumber triads where 
the smallest wave-number is greater than one-half of the middle wavenumber. 

In this subsection we will compute the shell-to-shell energy transfer in turbulence using field-theoretic method [^3. 
The procedure is identical to the one described for energy fluxes. Recall that the energy transfer rates from the m-th 
shell to the n-th shell is 

T nm =J2J2 ^"(k'lplq). (101) 

The p-sum is over m-th shell, and the k'-sum is over n-th shell fSection lTTTjl . The terms of S"s are the same as in flux 
calculation, however, the limits of the integrals are different. The shells are binned log arithmically with n-th shell 
being (fcos™ -1 , koS n ). We nondimensionalize the equations using the transformation 30] 

a a a 

k — — ; p = —v] q = —w, (102) 

u u u 

where a — kg-s"^ 1 . The resulting equation is 

,3/2 4Srf-l 



r 1 du r is "* 


- + 1 rl+v 




dv / 


/s _1 u Jus m - 


" J 1 V 



ii - K ^wws d L t L- *y M *-(~>"<*«>"*<».»>. v<*> 

whereF(v,w) was computed in the previous section. The renormalized parameters v* , and Kolmogorov's constant 
Kko required to compute T nm /U are taken from the previous calculations. From Eq. 110311 we can draw the following 
inferences: 



1. The shell-to-shell energy transfer rate is a function of n — m, that is, <I> nm = <&( n -i)(m-i)- Hence, the turbulent 
energy transfer rates in the inertial range are all self-similar. Of course, this is true only in the inertial range. 

2. T nn /n = 0. 

We compute the integral of Eq. IIO.'H and substitute the value of Kko, which yields T nm . The plots of T„ m for 2D 
and 3D fluids are shown in Fig. HJ In 3D the energy transfer is forward and local. In 2D however the energy transfer 
is forward for the nearest neighbours, but it is backward for fourth neighbour onward; these backward transfers are 
one of the major factors in the inverse cascade of energ y |5 fl. The sum of all these transfers is negative energy flux, 
consistent with the inverse cascade result of Kraichnan |24j. For details refer to Verma et al. [H^ |. 

Verma et al. computed the shell-to-shell energy transfer in 3D fluid turbulence using numerical simulations. 
Their result is shown in Fig. Comparison of Fig. Eland Fig. shows that theoretical and numerical computation 
of shell-to-shell energy transfer are consistent with each other. 

Incompressible fluid turbulence is nonlocal in real space due to incompressibility condition. Field-theoretic cal- 
culation also reveals that mode-to-mode transfer S(k\p\q) is large when p <§C k, but small for k ~ p ~ q, hence 
Navier-Stokes equation is nonlocal in Fourier space too. However, in 3D shell-to-shell energy transfer rate T nm is 
forward and most significant to the next-neighbouring shell. Hence, shell-to-shell energy transfer rate is local even 
though the interactions appear to be nonlocal in both real and Fourier space. Refer to Zhou [fi2l |. Domaradzki and 
Rogallo [lfj], Verma et al. yjfj, and Verma 

With this we conclude our discussion on shell-to-shell energy transfer in fluid transfer. 



C. EDQNM Calculation of Fluid Turbulence 



Eddy-damped quasi-normal Markovian (EDQNM) calculation of turbulence is very similar to the field-theoretic 
calculation of energy evolution. This scheme was first invented by Orszag [39] for Fluid turbulence. 
The Navier-Stokes equation is symbolically written as 

' p+q=k 
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Figure 6: Plots of shell-to-shell energy transfer rates T nm /H vs. n — m for 3D and 2D fluid turbulence. 
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Figure 7: Plot of normalized shell-to-shell energy transfer Tnm/n vs. n — m for d = 3obtained from numerical simulations on 
512 3 DNS. The nth shell is (k s n : fcos n+1 ) with s — 2 1 / 4 .The energy transfer is maximum for n = m ± 1, hence the energy 
transfer is local and self-similar. The energy transfer is also forward. Taken from Verma et al. [5^ |. 



where X stands for the field u, X(p,t)X(q,t) represents all the nonlinear terms, and C is the dissipation coefficient 
{y) . The evolution of second and third moment would be 



(1 + 20^ (x(M)x(-k,t)) 
{jt + C(P +p2 + q2) ) {X{ ~ K t)x{p ' t)x(q ' t]) 

If X were Gaussian, third-order moment would vanish. How 
correlation; here we replace (XXXX) by its Gaussian value, 



= W-M)X(p,t)X(q,t)) (104) 

p+q=k 

(X(q,t)X(p,t)X(v,t)X( S ,t)) 

p+q+r+s=0 

;ver, quasi-normal approximation gives nonzero triple 
which is a sum of products of second-order moments. 
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Hence, 

(X(-k,t)X(p,t)X(q,t)) = f dTexp-(;(k 2 +p 2 + q 2 )(t-T) V 

J ° P+q=k 

[<X (q, r)X(-q, r)) (X(p, r)X(p, r)> + ...], 

where ... refers to other products of second-order moments. The substitution of the above in Eq. H104|) yields a closed 
form equation for second-order correlation functions. Orszag '39] discovered that the solution of the above equation 
was plagued by problems like negative energy. To cure this problem, a suitable linear relaxation operator of the 
triple correlation (denoted by fx) was introduced (Eddy-damped approximation). In addition, it was assumed that the 

characteristic evolution time of (XX) (XX) is larger than (pkpq + K^ 2 + P 2 + Q 2 )) (Markovian approximation). As 
a result the following form of energy evolution equation is obtained 

(j t + 2(kA (X(k, t)X(-k, t)) = J dp9 kpq (t) J2 [<*(<!> *) (X(p, t)X(-p, t) + ...], (105) 

p+q=k 

where 



Okpq(t) = (1 - CXp -(fl k + [l v + fl q )t) I (/ife + [l v + fi q ) 

with 

fi k = [v + ri)k 2 + C s ^jf dq(E u {q))q^j . (106) 

The first and second terms represent viscous and nonlinear eddy-distortion rates respectively. Note that homogeneity 
and isotropy are assumed in EDQNM analysis too. 

The right-hand side of Eq. H105f) is very similar to the perturbative expansion of S uu (k\p\q) (under t — > oo). The 
term fi k of Eq. IllOfiH is nothing but the renormalized dissipative parameters. Thus, field-theoretic techniques for 
turbulence is quite similar to EDQNM calculation. There is a bit of difference however. In field-theory, we typically 
compute asymptotic energy fluxes in the inertial range. On the contrary, energy is numerically evolved in EDQNM 
calculations. 
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IX. INTERMITTENCY IN FLUID TURBULENCE 



The famous Kolmogorov's turbulence model assumes a constant energy flux or dissipation rate at all scales, i. e., 
H(k) ~ j dkv(k)k 2 E(k) is independent of A:. The renormalized viscosity v{k) ~ fc -4 / 3 and E(k) ~ fc -5 / 3 are consistent 
with the above assumption. Landau however [2H pointed out that the dissipation rate, which is proportional to the 
square of vorticity, is singular and quite inhomogeneous. Thus Kolmogorov's theory of turbulence needs modification. 
The above phenomena in which strong dissipation is localized both in time and space is called intermittency. 



A. Quantitative Measures of Intermittency 

There are several quantitative measures of intermittency. Consider the increment of the velocity, or some other 
field, between two points separated by 1, 

<5u(x,l) = u(x + l)-u(x). (107) 

The longitudinal component of <5u(x, 1) will be given by 

Su\\(l) =<Ju(x,l).l/Z, (108) 

and the transverse component is Su±(l) = <5u(x, 1) — 5un(l)l/l, Here we have assumed homogeneity and isotropy for 
turbulence, so that the increment in velocities depend only on I, not on x. Now we define longitudinal and transverse 
structure functions using 

S^(l) = ([S Ul] (l)] n ), U<ri(l) = ([6uj_(!)] n ) (109) 

respectively. The structure function S^{1) is expected to have a power law behaviour for I in the inertial range, 

SW{l)=aJ<", (110) 

where a n and (, n are universal numbers. The exponents £„ are called the intermittency exponent. 

Moments and probability density function (pdf) are equivalent description of random variables. Note that if 
P(8u\\(l)) were gaussian, i. e., 



P(6 U]l (l)) = — ^=exp- V ll J 2 " (111) 



1 (^n(O) 2 

— t= <■- 
ay v 71 " 

then, it is easy to verify that 

((6u\\(!)) n )<xo?. (112) 

Kolmogorov's model of turbulence predicts that 

or - e 1/3 / 1/3 , (113) 

For constant e, we obtain 

s w (i)«£ n/3 r /3 . (ii4) 

Systems with gaussian probability distribution or equivalently S^ n ' (I) oc l cn (c = constant) are called non-intermittent 
system. For intermittent systems, the tails of pdf decays slower that gaus sian, and could follow a powerlaw. 
Structure function can be written in terms of local dissipation rate [38] 

5ui - e 1 / 3 ! 1 / 3 . (115) 
Kolmogorov introduced the refined similarity hypothesis relating structure function to e; as 

S-f")(0=d n (er /3 )/ n/3 - (116) 

If 

(er>~l"», (117) 
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then 

Cn = |+/V3- ( 118 ) 

Many researchers have attempted to model e/. 

In any numerical simulation or experiment, the powerlaw range is quite limited. However, when we plot S^ n \l) 
vs. S( 3 >(1), we obtain a much larger scaling range. This phenomena is called Extended self-similarity (ESS). Since, 
SW(l) oc / (23, Cn measured using Eq. IjllOfl or ESS are expected to be the same. 

There have been some ingenious attempts to theoretically compute the intermittency exponents (e.g., see series 
of papers by L'vov and Procaccia (^3) . Yet, this problem is unsolved. Here we list some of the phenomenological 
models. For extensive discussions on intermittency refer to Frisch (T^ |. 



B. Results on Intermittency in Fluid Turbulence 

In fluid turbulence, the pdf of velocity increment deviates from gaussian In experiments and simulations one 
finds that Cn vs. n is a nonlinear function of n. Hence, fluid turbulence shows intermittency. Note that C2 ~ 0.71, 
which yields a correction of approximately 0.04 to Kolmogorov's spectral index of 5/3. However, the correction for 
large n is much more. See Frisch [lg for further details. 

Remarkably, starting from Navier-Stokes equation, Kolmogorov |2JJ| obtained an exact relation 

Sf| 3) (0 = -~d (119) 

under v — > limit (also see 0, H3|). Note that e is the mean dissipation rate. Unfortunately, similar relationship 
could not be derived for other structure functions. In the following discussion we will discuss some of the prominent 
intermittency models for fluid turbulence. 



1. Kolmogorov's log-normal model 

Obukhov and Kolmogorov jl^ claimed that the dissipation rate in turbulent fluid is log-normal. As a 
consequence, 

n n(n — 3) , 

where 

(e(x)e(x + l)>~r". (121) 

Numerical simulations and experiments give /x » 0.2. 

The predictions of this model agree well with the experimental results up to n « 10, but fails for higher values of 
n. In Fig. Owe have plotted the above Cn along with other model predictions given below. 



2. The (3- model 

Novikov and Stewart and Frisch et al. proposed that smaller scales in turbulent fluid is less space filling. 
In each step of the cascade an eddy Su n of scale l n splits into 2 D (3 eddies of scale l n+ \ = l n /2, where D is the space 
dimensionality, and (3 is a fixed parameter with < (3 < 1. In this model 

n S 

C»= 3 -3(^-3), (122) 

where /3 = 2~ s . 

Note that Cn is linear in n, and it does not match with experimental and numerical data for large n (see Fig. |HJ. 
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Figure 8: Plots of £ n vs. n for various intermittency models in fluids and MHD. She-Leveque's log-Poission model fits best 
with the experimental data in both fluid and MHD. For MHD turbulence, Kolmogorov-like models are in better agreement 
than KID's like model. 

3. The Multifractal Models 

Parisi and Frisch [4f| developed a multifractal model of turbulence. Maneveau and Sreenivasan [32j constructed an 
intuitive model. Here the energy cascade rate e/ is not distributed equally into smaller eddies, say, in each cascade 
it gets divided into pti and (1 — p)ei. After several cascades, one finds that energy distribution is very skewed or 
intermittent. The intermittency exponent in this model is 



with 



D„ 



C„ = (~ - l) D n + 1, 



log 2 ( P n + (l-p) n ) 1/{1 - n) 



(123) 
(124) 



For p near 0.7, ( n fits quite well with the experimental data. The deficiency of this model is that it requires an 
adjustable parameter p. For more detailed discussion, refer to Stolovitzky and Sreenivasan [Ifij. 



4- The Log-Poisson Model 



She and Leveque 43] proposed a model involving a hierarchy of fluctuating structures associated with the vortex 
filament. In their model 



Cn = !(l-£)+C (l-/?"/ 3 ) 



where Co is co-dimension of the dissipative eddies, and x and f3 are parameters connected by 

x 



C = 



1-/3 



(125) 



(126) 



For Kolmogorov scaling, x = /3 = 2/3. In hydrodynamic turbulence, the dissipative eddies are vortex filaments, i.e., 
one-dimensional structures. Therefore, the co-dimension is Co = 2. Hence, for fluid turbulence 



C SL = - + 2 

«m g f 



n/3 



(127) 
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The above prediction fits remarkably well with experimental results. All the above functions have been plotted in 
Fig. |H|for comparison. 

X. CONCLUSION 

In this paper we review some of the important results in statistical theory of fluid turbulence. The focus was on 
field-theoretic techniques applied to fluid turbulence. 

Many current topics of interest have not been included in the paper. Some of them are intermittency, large eddy 
simulations, compressible turbulence etc. 
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Appendix A: FOURIER SERIES VS. FOURIER TRANSFORM FOR TURBULENT FLOWS 

In statistical theory turbulence we typically assume the flow field to be homogeneous. Therefore, Fourier transform 
is not applicable to these flows in strict sense. However, we can define these quantities by taking limits carefully. This 
issue has been discussed by Batchelor Q and McComb jgj. We briefly discuss them here because they form the basis 
of the whole paper. 

A periodic function u(x) in box L d can be expanded using Fourier series as following: 



u(x) = V" u (k) exp (ik • x) , (Al) 

u(k) = Jd J dxu ( x ) CX P (-* k ' x ) > ( A2 ) 

where d is the space dimensionality. When we take the limit L — > oo, we obtain Fourier transform. Using u(k) = 
u(k)L d , it can be easily shown that 

f dk 

U(X) = J (2^p u(k)exp(?k,x) ' (A3) 
u(k) = J dxu (x) exp (-ik • x) , (A4) 

with integrals performed over the whole space. Note however that Fourier transform (integral converges) makes sense 
when u(x) vanishes as |x| — > oo, which is not the case for homogeneous flows. However, correlations defined below 
are sensible quantities. Using the above equations, we find that 



(ttj(k)uj(k')) = / dxdx' (ui(x)Uj (x')) exp — i(k-x + k' • x') 

= J drCij(r) exp — ik • r J dxexp — z(k+k') • x 

= C lJ (k){2n) d 5{k + k r ) (A5) 

We have used the fact that S(k) sa L d / (2ir) d . The above equation holds the key. In experiments we measure correlation 
function C(r) which is finite and decays with increasing r, hence spectra C(k) is well defined. Now energy spectrum 
as well as total energy can be written in terms of C (k) as the following: 



<« 2 > = ^ / ^ 2 = E |G(k)| 2 -^j^y (Kk)f 

= ( ^ 1} /(^ C(k) (A6) 



We have used the fact that 5(k) w L d /(2ir) d . Note that ^|u(k)r^ = (d - l)C{k)L d [see Eq. (SB] is not well defined 
in the limit L — > oo. 

In conclusion, the measurable quantity in homogeneous turbulence is the correlation function, which is finite and 
decays for large r. Therefore, energy spectra etc. are well defined objects in terms of Fourier transforms of correlation 
functions. 

We choose a finite box, typically (27r) d , in spectral simulations for fluid flows. For these problems we express the 
equations (incompressible MHD) in terms of Fourier series. We write them below for reference. 

fj^ + vk 2 ^j u*(k,t) = -ikjptot (k, t) - ikj y^[Mj(q, t)uj(p, t)] (A7) 

The energy spectrum can be computed using Uj(k, t): 

J Emk = £ /,_/*, ,2 - / dk kwi> ft m 
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where n is the lattice vector in (i-dimensional space. The above equation implies that 

E(k) = \^l Sd k d -\ (A9) 

A natural question is why the results of numerical simulations or experiments done in a finite volume should 
match with those obtained for infinite volume. The answer is straight forward. When we go from size 2ir to L, the 
wavenumbers should be scaled by (2ir)/L. The velocity and frequency should be should be scaled by (2ir)/L and 
[(2n)/L} 2 to keep dimensionless v fixed. The evolution of the two systems will be identical apart from the above 
factors. Hence, numerical simulations in a box of size 27r can capture the behaviour of a system with L — > oo, for 
which Fourier transform in defined. 



Appendix B: PERTURB ATIVE CALCULATION OF NAVIER STOKES EQUATION 

The Navier Stokes can be written as 

Ui(k) = G(k) - l -P+ m (k) J dp[ Uj (p)u m (k p)] (Bl) 

where the Greens function G can be written as 

G-\k,Lj) = -iuj-Z uu (B2) 

We solve the above equation perturbatively keeping the terms upto the first nonvanishing order. The corresponding 
Feynmann diagram is 




(B3) 



The solid line represents fields u, and the wiggly line (photon) denotes G. The filled circle denotes — (i/2)P^ 
vertex. These diagrams appear in renormalization calculations as well as in energy flux calculations. 

1. Viscosity Renormalization 

The expansion of u in terms of Feynman diagrams are given below: 




Factor of 2 appears in I u because of <> symmetry in the corresponding term. To zeroth order, the terms with <> 
are zero because of quasi-gaussian nature of > modes. To the next order in perturbation, the third term of I u is 
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o-+ o- 

fV. 



- + 



+higer oder diagrams (B5) 

In the above diagrams solid lines denote < u(k)ii(k') >. As mentioned earlier, the wiggly line denotes Green's 
functions. All the diagrams except 4,8,12, and 16th can be shown to be trivially zero using Eqs. 1731 7411 We assume 
that 4,8,12, and 16th diagrams are also zero, as usually done in turbulence RG calculations [23, Ltii IlLt L2j| ■ Hence, 
the term is zero. Now we are left with >> terms (3rd term of I u ), which is 



q = < =-SX(k) (B6) 



where 



(d-l)*E = « + (B7) 



o 



In the above equation we have omitted all the vanishing diagrams (similar to those appearing in Eq. |B5| ). These 
terms contribute to Eg. 

The algebraic expressions for the above diagrams are given in Section IVlTl For isotropic flows, the algebraic factors 
S{k,p,q) resulting from tensor contractions are given below. 

S(k,p, q) = P b + jm (k)P+ ab (p)P ]a (q) = kp ((d - 3)2 + 2z 3 + (d - l)xy) (B8) 

In the next subsection we will derive the terms for mode-to-mode energy transfer function. 

2. Mode-to-Mode Energy Transfer in fluid Turbulence 

In Section 3, we studied the "mode-to-mode" energy transfer S , ""(k;'|p|q) from mode p to mode k', with mode q 
acting as a mediator. The perturbative calculation of S involves the following terms 

(S uu (k'\p\q)) = \ + ^^^^ + t # (B9) 

In all the diagrams, the left vertex denotes fcj, while the filled circle of the right vertex represent (—i/2)Pjj m - For 
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isotropic nonhelical flows, the algebraic factors are given below. The factors for the diagrams are Ti, T2, T 3 in sequential 
order. 

Ti(A,j»,?) = k l P+ b (k)P ja {p)P lb {q) = kp((d-3)z + (d-2)xy + 2z 3 + 2xyz 2 + x 2 z) (BIO) 
T 3 (k,p,q) = -k t P+ b (p)P ia (k)P lb (q)=-kp({d-3)z + {d-2)xy + 2z 3 + 2xyz 2 + y 2 z) (Bll) 
T 3 (k,p,q) = ~k l P+ b (q)P ja (k)P ]b (p)=-kq(xz-2xy 2 z-yz 2 ) (B12) 

(B13) 

These terms are similar to those given in Leslie [3JJ|. 

Appendix C: MODE-TO-MODE ENERGY TRANSFER IN SCALAR TURBULENCE 

The equations for passive scalar turbulence are 

— + (u-V)u = -V;p + i/V 2 u, (CI) 

— + (u-V)0 = kV 2 0, (C2) 
Vu = 0, (C3) 

where 9 is the scalar density. Note that the scalar is convected by velocity field, but the scalar does not affect the 
flow. For details refer to Lesieur J^U and Stanisic [4fj |. 

For energy transfer in scalar field, we can follow the same procedure as in fluid turbulence. If we take only a single 
triad (k 7 . p, q) with k' + p + q = 0, energy is conserved for k = 0. Following fluid turbulence, we can show that the 
energy equation for scalar turbulence can be written in terms of 'mode-to-mode energy transfer' 5 88 (k'|p|q) from 
mode 0(p) to 0(k') with 0(q),u(q) as a mediator, which is 

5"(k'|p|q) = -9([kU(q)][0(fc')0(p)]) (C4) 
where 9 stands for the imaginary part of the argument. The energy equation for scalar field is 

d 



j ft +2Kk 2 )C e (k) = [S e9 (k'|p|q)+<^(k'|q|p)] (C5) 

Note that there is no cross-transfer between u and i\> energy. It is also important to note that both C u and are 
conserved in every triad interaction, i.e., 

S«"(k'|p|q) + 5™ l (k'|q|p) + S" u (p|k'|q) + 5™ l (p|q|k') + S""(q|k'|p) + S""(q|p|k') = (C6) 

S^(kqp|q) + ^(k / |q|p) + 5^(p|k'|q) + 5^(p|q|k')+^(q|k'|p)+5^(q|p|k') = (C7) 

These are the statements of "detailed conservation of energy" in triad interaction (when v = K = 0) 
The energy flux 11^ from a wavenumber sphere of radius ko is Q 

^-LhU 1 ^ 1 <C8) 



Appendix D: MODE-TO-MODE ENERGY TRANSFERS IN RAYLEIGH BENARD CONVECTION 

One of the most studied model of convection is Rayleigh Benard (RB) Convection. In this model the fluid con- 
fined between two parallel plates is heated from below. Constant temperature is maintained across these plates. 
Nondimensionalized equations for RB model are 

-^ + (u-V)u = -Vcr + RPOi + PV 2 u, (Dl) 
— + (u-V)9 = u 3 + V 2 6, (D2) 
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where the nondimensinal parameters R and P — v/k are Rayleigh and Prandtl numbers respectively. The fluid is 
assumed to incompressible except for the buoyancy term. 

We can apply 'mode-to-mode energy transfer' model to RB convection. The above hydrodynamical equations in 
Fourier space are 



ihcrk-ikj u M) u i(p) +RP9(k)6ia -Pk 2 u t {k) (D3) 



du t {k) 

at 

d9(k) 

at 

k lUl {k) = 0, (D5) 



p+q=k 

u a (k)-ikj U MMP) - k 2 0(k), (D4) 

p+q=k 



where k = p + q. We can derive interesting results by focussing on a single triad (k', p, q) such that k' + p + q = 0. 
Clearly k' = k. We can easily derive the following energy equations: 



d l u (k)l _ uurtj\„\„\ , cmu/i./i_i_\ , r>T><xsmr*.\.. fi./M o nu 2 l u (k) 12 

at 2 



5 uu (k'|p|q) + 5 uu (k'|q|p) + RPn [9(k)u 3 (k')} - 2Pk 2 ^-^, (D6) 



d \0(k)l = s ee {k , lp]ci) + ^ (k '| q | p) + g [e(k)u 3 (k')} - 2fc 2 M^ (D7) 



where 



^"(k'|p|q) = -3([k'-u(q)][u(k')-u(q)]), (D8) 
5 ee (k'|p|q) = -3 ([k' • u(q)] [0(k') ■ 0(q)]) . (D9) 

Here and 3 represents the real and imaginary part of the argument. The quantity 5' ee (k'|p|q) represents the 
energy transfer from mode 6>(p) (the field variable with the second argument) to mode 0(k') (the field variable with 
the first argument) with the help of the mode 0(q) (the field variable with the third argument) acting as a mediator. 
The above formalism is called "mode-to-mode" formalism for energy transfer an d it differs from Kraichnan's 

" combined energy transfer" formalism |2.'il29l|. The energy equation can be interpreted as follows: The field variables 
with wavenumber k' [u(k'), #(k')] receives energy from the modes p and q through mode-to-mode energy transfer 
terms, and it also receives energy due to interaction term 0(k)tt3(k'). 

It is important to keep in mind that |6*(k)| 2 /2 is not real energy. However, it has structure of energy, and we can 
apply the energy transfer formalism here as well. This idea has been exploited heavily in the past to infer the direction 
of flux of passive scalar etc. [3]. Dar et al. Q and Verma have shown that the sum of all energy transfer rates 
along u-u and 9 — 9 channels are zero, i.e., 

S XX (k'|p|q) + S X *(k'|q|p) + S XX (p\k'\q) (D10) 
+S XX (p\ q \k) + S XX (q\k / \p) + S XX (<l\p\k / ) = 0, (Dll) 

where XX could be uu or 99. 

Using this identity we can easily show that without viscous and thermal diffusion 

^|u(k)| 2 + |u(p)| 2 + |u(q)| 2 = 2flPR[0(k)u$(k) + 0(p)u$(p) + 0(q)u$(q)]. (D12) 

| [|0(k)| 2 + |0( P )| 2 + |0(q)| 2 ] = 25R[0(k)u*(k) + 0(p)«5(p) + tf(q)«S(q)], (D13) 

The interpretation of the above equations is that the triads [u(k'), u(p), u(q)] and [6*(k'), 9(p), 9(q)) exchange energy 
between each other via 9(k)u 3 '(k) interaction terms. The mode-to-mode interaction conserves energy within a triad. 
The viscous and diffusive terms dissipate kinetic energy and ^-energy respectively. 
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